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Abstract 



An analysis of local fluctuations, or spikes, is performed for charged particles 
produced in central C-Cu collisions at 4.5 GeV/c/nucleon. The distributions 
of spike-centers and the maximum density distributions are investigated for 
different narrow pseudorapidity windows to search for multiparticle dynami- 
cal correlations. Two peaks over statistical background are observed in the 
spike-center distributions with the structure similar to that expected from the 
coherent gluon radiation model and recently found in hadronic interactions. 
The dynamical contribution to maximum density fluctuations are obtained 
to be hidden by statistical correlations, though behavior of the distributions 
shows qualitative agreement with that from the one- dimensional intermittency 
model. The observed features of the two different approaches, coherent vs. 
stochastic, to the formation of the local dynamical fluctuations are discussed. 
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1 Introduction 



Reflecting an intermittent structure of distributions of particles produced in 
high-energy collisions, local fluctuations play an important role in hadroniza- 
tion process. Intermittency/fractality studies |IJ show these fluctuations to 
be significant in investigations of multihadronic production as well as to be a 
possible signal of quark-gluon plasma formation. Although a vast activity in 
study self- similarity have disentangled contributions of different known mech- 
anisms to production of dynamical fluctuations/correlations, an origin of the 
latters remains still unclear. This arises need for applying more direct (and 
additional) methods of studying properties of dense groups of particles. 

The data used here has been considered in our recent studies of self-similar 
nature of multihadronic production ||. Strong multifractality has been 
found, pointing out a possible non-thermal phase transition and two different 
regimes in particle production during the cascade. This conclusion confirms 
results of our earlier studies [§], f| made with less statistics and has been re- 
cently supported by similar observations reported for ultra-relativistic nuclear 
collisions . Further analysis || have shown more structure in dynamical fluc- 
tuations indicating chaotic nature of multiparticle production with a specific 
scaling-law referred to as erraticity J7|. 

Although the self-similar underlying dynamics has been obtained in searches 
in fractality terms, no dynamical effects was revealed in the investigations of 
maximum density distributions ||. However, the analogous study performed 
for hadron-emulsion data in the energy range of 200-400 GeV || has indicated 
the existence of dynamical multiparticle correlations and clustering seen for 
high-density spikes. 

In this letter we investigate local pseudorapidity fluctuations, or spikes, an- 
alyzing two types of distributions, assigned to coherent vs. chaotic approaches 
in particle-emission process. Respectively, the distributions of centers of spikes 
and maximum density distributions are considered, based on the methods, 
used earlier and shown to be distinguished by the reliability and the high 
certainty in searching for multiparticle dynamical correlations JiTi|-|i"5||. Study 
of fluctuations in limited regions of collision phase-space allows reducing of 
contribution from conservation constraints and elicits underlying dynamics of 
particle-production process. 

2 Data sample and analysis procedure 
2.1 Data sample 

The results presented are based on experimental data came from interactions 
of the JINR Synchrophasotron (Dubna) 4.5 A GeV/c 12 C beam with a copper 
target inside the 2m Streamer Chamber SKM-200 ||14|| . A central collision 
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trigger was used: absence of charged particles with momenta p > 3 GeV/c in 
a forward cone of 2.4° was required. A more detailed description of the set-up 
design and data reduction procedure are given elsewhere fl4] , |l5|l . Systematic 
errors related to the trigger effects, low-energy pion and proton detection, the 



admixture of electrons etc. does not exceed 3% [|T5| .|T6 



The scanning and the handling of the film data were carried out on special 
scanning tables of the Lebedev Physical Institute (Moscow) ]17|]. The average 
measurement error in the momentum (e p /p) was about 12%, and that in the 
polar angle measurements was (e&) ~ 2°. In total, 663 events with charged 
particles in the pseudorapidity range of Ar] = 0.2 — 2.8 (in the laboratory 
frame) were processed. The angular measurement accuracy does not exceed 
0.1 in the 77-units. In addition, particles with p^ > 1 GeV/c are excluded from 
the investigation as far as no negative charged particles were observed with 
such a transverse momentum. Under the assumption of an equal number of 
positive and negative pions, this cut was applied to eliminate the contribution 
of protons. After the kinematic cuts, the mean multiplicity is 23.8 ± 0.4. 



2.2 Analysis procedure 

The multiparticle fluctuations are studied in the pseudorapidity phase-space 
regions. The following procedure is used to search for dynamics of the fluc- 
tuations. For each event the ordered pseudorapidities, rj = — In tan (■$ is 
the polar angle of the particle), are scanned with a fixed pseudorapidity win- 
dow (bin) across the full //-range of the event, and the spike with 5n number 
of tracks, hit in the window 5r], is determined. Then the centers of spikes, 
rjo = (l/5n) J2^=iVji are calculated for all events and the distribution in r] is 
investigated to reveal dynamical correlations. The distribution with respect to 
the maximum density fluctuations, defined as p max = 6n max /5r], where Sn max 
is the maximum number of particles per spike in each event for the chosen Sr), 
is analyzed as well in order to obtain dynamical character of the fluctuations 
observed. 

The conclusions about dynamical content of the fluctuations could be af- 
fected by the dependence of such an analysis on the form of pseudorapidity 
distribution and the fact that the method of maximum fluctuations deals with 
the mostly populated regions of the distribution. To avoid these and to com- 
pare the results from different experiments, the "cumulative" variable, 

77(77) = / p( V ')dr,'/ / p(r/)dr/, (1) 

with the uniform spectrum p(v) within the interval [0,1] was proposed to be 
used |TJ|. The transformed variable ([]]) are usually considered in studying in- 
termittency |J, also utilized in our recent reports Note, pseudorapidity 
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is argued to be the most suitable variable to analyze correlations providing in- 
termittent structure of high-energy events |H| . 



3 Results and discussion 
3.1 Spike-center distributions 

The pseudorapidity spike-center distributions for the different size 5 ?7-bins and 
for spikes of the different density are presented in Fig. 1. The widths of the 
bins are chosen to be compatible with those used earlier P, p0| . Multi-peak 

structure of the distributions can be seen for 5 V= 0.04 (Srj ~ 0.1) (Fig. la) 
and 0.08 (~ 0.2) (Fig. lb). However, one can observe two peaks placed in 
the region about the same ?7crP os iti ons with a tendency of the distributions 
to have a double-peak shape as the size of the bin increases. The two peaks 
become much more pronounced when the events are scanned by the large 5 V, 
e.g. of the width of 0.2 as shown in Fig. Id. Fitting these two bumps with 
Gaussians, the peaks averaged over the different spikes, are found to be placed 
at 0.17 and 0.57. Recounted to the i]- variable, the positions of the peaks are 
of the values of 0.60±0.05(stat) ±0.12(syst) and 1.30±0.03(stat) ±0.10(syst) 
with the distance, 

d = 0.68 ± 0.06(stat) ± 0.16(syst) (2) 

between them. 

Recently, the study of the spike-centers has been carried out for hadronic in- 
teractions at 205-306 GeV/c pi] . The double-peak shape of the ?7o-distribution 
for pp-collisions vs. a single-peak structure in 7r/Kp-interactions have been ob- 
served in agreement with the coherent gluon-jet emission model [|H]]. The peaks 
has been found to be separated by the distance of 0.57±0.03(stat) ±0.12(syst), 
also consistent with the model prediction. The double-peak form obtained here 
for the central nuclear interactions are similar to that for pp-type reaction, in- 
dicating superposition of nucleon-nucleon interactions in nucleus-nucleus one. 
Moreover, the value of the distance found exceeds that in hadronic collisions 



being in agreement with theoretical expectations |22[ . 

To observe dynamical correlation effect in these distributions, the analo- 
gous distributions have been obtained from the simulated pseudorapidity one- 
particle spectrum p( ? 7), in which, evidently, any information of two or more 
multiparticle correlations is lost. The simulation procedure was as follows. In 
accordance with the multiplicity distribution of the experimental sample we 
have randomly generated corresponding number of tracks. Then the pseudora- 
pidities has been distributed according to the real p(??)-spectrum in a quantity 
of the generated multiplicity. The total number of the events simulated was 
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66300, so that it exceeded the experimental statistics by a factor of 100. It is 
clear that the statistical properties of this set are completely analogous to those 
of the ensemble resulted from arbitrary mixing of tracks from different events, 
subject to the condition of retention of the p(^)-distribution, and the obtained 
sample represents the result of independent particle emission hypothesis. 

The ^-distributions of the simulated events are shown in Fig. 1 by open 
circles. Remarkable difference one can observe between the experimental dis- 
tributions (solid circles) and those obtained in assumption of completely mi- 
correlated particle production. No whatsoever peaks are seen in the latter case, 
following the background level and manifesting the double-peak structure to 
be the dominant one. 

From the comparison of the experimental and simulated distributions in the 
spike-centers proceeded for various 5 ^-windows and different many-particle 
spikes, one can apparently conclude about dynamical character of the produc- 
tion of high-density fluctuations. The dynamics of the intermittent structure 
of the data studied could be assigned to the model of coherent gluon radia- 



tion [10] as to one of a real candidate of the mechanism of formation of the 
fluctuations. 

To assess the reliability of the above conclusions and the stability of the 
results we vary the investigation with changing the Arrange under consider- 
ation as well as the polar angle, within the experimental error (£■&). The 
observed character of the distributions remained unchanged. 

3.2 Maximum density fluctuations 

The further analysis deals with the normalized distributions, (l/N)dN/dp max , 
of iV events as shown in Fig. 2 for four different scanning windows. 

The exponential decrease at p max > (p max ) of these distributions, averaged 
over all multiplicities n, is analogous to that observed in our previous stud- 
ies H^(J and in other investigations performed for different reactions |9|. p3 



Such a behavior differs from the poissonian one, expected for processes with 
weak correlations of produced hadrons or for models of multiperipheral or 
Regge types taking into account a limited number of reggeons. The exponen- 
tial behavior of the p max -spectra is argued to be a consequence of primordial 
multiparticle correlations which are irreducible to two-particle correlations |13| . 
The given inequality between the dispersion and the mean values (p ma x) con- 
firms the non-poissonian character of the distributions, pointing out significant 
contribution of multi-particle correlations to the local fluctuations observed. 

Although the transition to the "cumulative" variable did not influence much 
the above conclusion compared to the earlier results |8],[21J, it seems to be 



essential in studying the distributions at high p max -values. Indeed, if the bell- 
like shape for small bins (Fig. 2a) is similar to that observed in the previous 
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reports, an increase of the 5 ?7-width (Fig. 2b, c) makes large p max -tails to 
appear in the p max -spectra. 

A change of the shape of the maximum density distribution from the ex- 
ponential to more flat one with increase of the p max seems to be in agreement 
with that expected from the one-dimensional intermittency model ||11|| . The 
key feature of the model is an existence of two regimes in particle produc- 
tion process, described as turbulent and laminar components, leading to two 
maxima in the p max -distributions. 

It is worthwhile that the model makes predictions for the p max -distributions 
considered at given multiplicity n. Note, the fixed-n distributions are energy 
and reaction-type independent that allows increasing statistics by compiling 
results from different experiments. A study of the fixed-n distributions has 
been carried out for hadronic interactions and the large-p max tails have been 
observed for St] = 0.1 fl24 . 

In Fig. 2 we present the p max -distributions for two different narrow n- 
intervals (solid squares and triangles). An effect of deviation of the p max - 
distribution from the exponential behavior at large maximal densities are well 
seen and became more prominent compared to the case of mixed-n distribu- 
tions. Already at 5r) ~ 0.1 (Fig. 2a) the shape of the distributions develops 
tails at p max > (p max ), apparent also for larger 5 "-bins, 5 V= 0.12 and 0.2, 
but only for high-multiplicity events, 24 < n < 30 (triangles in Fig. 2b, c). 

For the widest 5 V shown, 5 V— 0.4 (Fig. 2d), one can not see the dif- 
ference in the behavior of the distributions for fixed n-regions. However, the 
distribution for all multiplicities shows slight flattening at large p max , but, in 
our opinion, this is rather due to the averaging over all n's available. 

Similar to the mixed-n spectra the fixed-n interval distributions show their 
non-poissonian character resulting to an inequality between the dispersion and 
the (pmax) and indicating contribution of multiparticle correlations. 

To reveal dynamical correlation effect, the obtained distributions are com- 
pared to those based on the above described sample of the simulated events, in 
which no whatsoever dynamical correlations exist. The resulted distributions 
for the four bins studied are shown in Fig. 2 for total multiplicity as well 
as for the two fixed-n intervals. The values of x 2 /DOF give a quite good fit 
between the data and the generated spectra independent of the width of the 
bins or the multiplicity. From this comparison one can not conclude about 
(intermittency) dynamics behind observed multiparticle fluctuations. It seems 
that the dynamical correlations in the distributions of maximum densities are 
too suppressed by statistical "noise" to be detected. 

However, it is worth to notice that the p max -distributions at 5 V— 0.12 and 
0.2 and 24 < n < 30 (fig. 2b, c), obtained from the uncorrelated emission 
hypothesis deviate from the data-based spectra at high-p max values: if the 
flattening is well-seen for the measured distributions, there is no change in the 
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behavior of these distributions from the simulated sample. 

The non-poissonian behavior found to be stronger for smaller bin-sizes and 
the deviations between the data and the simulated distributions observed for 
higher maximum densities, both seem to be a reason of non-statistical fluctua- 
tions found by us via the method of normalized scaled factorial moments [Q-|5| 
and leading to extreme fluctuations recently proposed to be searched for . 

All these indicates a need of more detailed analysis to be done with higher 
multiplicity reactions. Note, that such an analysis could be compared with 
that performed here due to the "cumulative" variable used and fixed-n intervals 
considered. 

A study of influence of the error (e&) in the measurement of the polar 
angle $ of the produced charge particles demonstrated stability of the obtained 
distributions and, therefore, the reliability of the conclusions above. 



4 Conclusions 

In summary, a study of spike production in central C-Cu collisions at 4.5 
GeV/c per nucleon is presented. To avoid the dependence of the results on 
the form of the pseudorapidity distribution the transformation to the uniform 
spectrum variable is utilized. The spike-center distributions and those in max- 
imum density fluctuations are investigated for various narrow pseudorapidity 
windows to obtain dynamical collective effects. 

The double-peak shape of the spike-center distributions is found, similar to 
that recently obtained in high-energy pp-interactions. The distance between 
the positions of the peaks seems to be in agreement with the expectation of the 
coherent gluon jet-emission model. Comparison with the results of completely 
uncorrelated particle-production model confirms dynamical origin of the effects 
observed. 

The maximum density distributions are investigated for mixed-n values as 
well as for fixed-n intervals. Non-poissonian character of the distributions is 
found indicating contribution of multiparticle correlations to the spikes. A 
flattening of the shape of the distributions at large maximum densities is ob- 
tained for high multiplicities in qualitative agreement with one-dimensional 
intermittency model. The dynamics are found to be hidden by strong statis- 
tical correlations, though visible deviations in the behavior of the simulated 
distributions compared to the measured ones are indicated. 

To conclude, a direct study of two different approaches - coherent vs. 
stochastic - to the formation of the local dynamical fluctuations in multiparti- 
cle production in central nuclear collisions at intermediate energy is performed. 
The coherent nature of dense group production is found to be clearly mani- 
fested with a structure similar to that recently found in hadronic interactions. 
Though no dynamics is obtained by means of the stochastic (intermittent) ap- 
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proach, the features of the maximum density distributions indicate a possible 
origin of self-similarity of the local fluctuations. 
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Figure captions 



Fig. 1. Experimental (•) and simulated (o) spike-center distributions for dif- 
ferent 5 ^-bins and various <5n-number of particles per bin: 

(a) 5 rj = 0.04, 5n = 4, (b) 5 rj = 0.08, 5n = 5, 

(c) S rj = 0.12, Sn = 7, (d) S rj = 0.2, 5n = 9. 



Fig. 2. Normalized experimental (solid symbols) and simulated (open sym- 
bols) distributions in maximum density p max for different 5 ^-windows and 
three multiplicity patterns: 

(a) 5 rj= 0.04, x 2 /DOF ~ 1.3 (all n, (p max ) ^ 39.4, RMS ~ 11.2), 
X 2 /DOF ~ 0.5 (14 < n < 20, (p max ) ~ 38.6, RMS ~ 8.1), 
X 2 /DOF ~ 0.7 (24 < n < 30, (p max ) ~ 49.2, RMS ~ 9.0), 

(b) 5 V= 0.12, x 2 /DOF ~ 1.2 (all n, (p max ) ~ 22.0, RMS ~ 6.6), 
X 2 /DOF ~ 1.3 (14 < n < 20, (p max ) ~ 12.8, RMS ~ 3.6), 
X 2 /DOF ~ 2.0 (24 < n < 30, (p max ) ~ 25.2, RMS ~ 4.4), 

(c) 5 rj= 0.2, x 2 /DOF ~ 1.1 (all n, (p max ) ~ 17.4, RMS ~ 5.2), 
X 2 /DOF ~ 1.0 (14 < n < 20, (p max ) ~ 14.0, RMS ~ 2.6), 
X 2 /DOF ~ 1.7 (24 < n < 30, (p max ) ~ 20.2, RMS ~ 2.8), 

(d) 5 V = 0.4, x 2 /DOF ~ 0.9 (all n, (p max ) ~ 13.2, RMS ~ 4.2). 
X 2 /DOF ~ 0.7 (14 < n < 20, (p max ) ~ 10.3, RMS ~ 1.6), 
X 2 /DOF ~ 0.9 (24 < n < 30, (p max ) ~ 18.2, RMS ~ 1.8), 
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Figure 2 
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